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ABSTRACT 

Galaxy source counts that simultaneously fit the deep mid-infrared surveys at 24 mi- 
crons and 15 microns made by the Spitzer Space Telescope and the Infrared Space 
Observatory (ISO) respectively are presented for two phenomenological models. The 
models are based on starburst and luminous infrared galaxy dominated populations. 
Both models produce excellent fits to the counts in both wavebands and provide an 
explanation for the high redshift population seen in the longer Spitzer 24 micron 
band supporting the hypothesis that they are luminous-ultraluminous infrared galax- 
ies at z=2-3, being the mid-infrared counterparts to the sub-mm galaxy population. 
The source counts are characterized by strong evolution to redshift unity, followed 
by less drastic evolution to higher redshift. The number-redshift distributions in both 
wavebands are well explained by the effect of the many mid-infrared features passing 
through the observation windows. The sharp upturn at around a milliJansky in the 
15 micron counts in particular depends critically on the distribution of mid-infrared 
features around 12 microns, in the assumed spectral energy distribution. 

Key words: Cosmology: source counts - Infrared: source counts. Surveys - Galaxies: 
evolution. 



1 INTRODUCTION 

The detection of the peak in the far-infrared background 
, at 140/im by the COBE-DIRBE instrument has empha- 
sized that much of the star formation in the Universe has 
taken place in dust enshrouded environments and that the 
galaxies responsible have undergone strong evolution to 
the present epoch (e.g. Hauser et al. (1998), Lagache & 
Puget 1(2000 )). This evolution has also been inferred and 
confirmed at wavelengths from 7-170/im by a progression 
of increasingly inspiring space missions from IRAS, ISO to 
Spitzer (Soifer et al. 119871 . Kessler et al. Ill99(ill . Werner et 
al. 42004) ) and at sub-mm wavelengths by the SCUBA in- 
strument IHolfand et af. 1999t . The emission in the dust 
spectrum of galaxies due to star formation peaks in the 
wavelength region between 60-100/im however a significant 
fraction (~40%) of a typical star-forming galaxy's luminos- 
ity is radiated at mid-infrared wavelengths from 7-40/im 
l |Soifer, Houck fc Neugebauer 1987| |. Moreover, unlike the 
far-infrared, observations at mid-infrared wavelengths are 
not so severely constrained by detector array size, resolu- 
tion and detector effects such as transients and glitches. The 
mid-infrared spectra of star-forming galaxies are dominated 
by a complex set of emission features at 3.3, 6.2, 7.7, 8.6, 
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11.3 & 12.7^m attributed to Polycyclic Aromatic Hydro- 
carbon (PAH) molecules and the silicate absorption feature 
at 9.7//m (Boulanger et al. (I1998L Peeters et al. J2004II ). 
These PAH features can contribute 10-30% of the bolomet- 
ric infrared luminosity of infra-red galaxies and are good 
indicators of the star- formation rate (Sp oon et al. 2004a| |. 
Typical widths of these features (up tolO^m) can be compa- 
rable to instrument filter bandwidths (e.g., 3.5/xm & 6/im for 
the ISOCAM-LW2(5-8.5Mm) filter and ISOCAM-LW3(12- 
18/im) filters respectively) and will therefore have a pro- 
found effect on any observations in the mid-infrared. 

The strongest constraint on the evolution in the mid- 
infrared until recently, has been the source counts at 15/im 
in the ISOCAM-LW3 band. The ISO counts sp an over 4 
orders of magnitude in fiux (lElbaz et af. 1999ll . The fiux 
density range from 150-0. 5mJy is covered by the European 
Large Area ISO Survey over four main fields over a to- 
tal of ~12sq.deg. (ELAIS, Oliver et al. (20O). At fiux 
densities >2mJy the source counts are Euclidean in na- 
ture and are well fitted by a combination of star-forming 
and quiescent galaxies with little or no contribution from 
more extreme luminous (LIG, Lir > 10^^ Lq) or ultralu- 
minous (ULIG, Lm > IO^^Lq) infrared galaxies (Serjeant 
et al. 1200011 . Gruppioni et al. 120021 Rowan- Robinson 
et al. 120041 '). At fluxes <2mJy the source counts show 
a strong departure from non-evolving scenarios. The fainter 
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ISO surveys detected extreme evolution producing super Eu- 
clidean slopes at fluxes < ImJy, flattening at fainter fluxes 
of 0.4m Jy (Fadda et al. ^M^ . Pozzi et al. lHHJ, Metcalfe 
et al. lUTH^ . Gruppioni et al. f2002j , Altieri et al. iTM^ . 
Aussel et al.. J1999t ') The majority of the deep ISOCAM 
sources have been identified with luminous infra-red galax- 
ies with luminosities of the order of 10^^ Lq (Flores et al. 
111999 ')'). The ISO 15/im counts are well fit by models as- 
suming a rapid increase in the star formation activity be- 
tween 0<z<l (e.g. Pearson J2UUlaL Xu et al. J2UU3L La- 
gache et al. J2003II '). The MIPS 24/im band on the Spitzer 
Space Telescope is well suited to the study of distant star 
forming galaxies at z>l since the 24/im band is sampling 
the rest frame PAH emission from 6-12/im from 0.7<z<2.5. 
Observations in the 24^m band over the flux range 50- 
0.06mJy have confirmed the strong evolution of the mid- 
infrared population (Papovich et al. (2111) I j. Marleau et 
al. (2004)), but however cannot be explained by evolution- 
ary models based solely on the previous ISO 15/im popula- 
tion (e.g. Xu et al. lEPOS), Chary & Elbaz (2003), King 
& Rowan-Robinson ll^IMIl . Lagache et al. ll200^ i. The 
ISO 15/im differential counts peak at around 0.4mJy with 
a median redshift ~0.8 (Elbaz et al. 11999II . Elbaz et al. 
1200211 1. conversely the peak in the 24pm counts occurs at 
fiux densities almost twice as faint around 0.2mJy imply- 
ing a higher redshift distribution for the dominant sources 
jChary et al. 2004j for a realistic range of the 24pm/15pm 
fiux ratio ^Gruppioni et al. 2005^ . 

In this letter the emphasis is on the simultaneous fit- 
ting of the source counts in the /S'0-ISOCAM-LW3 15pm 
and Spitzer-MIPS 24/im bands. Detailed results on the mod- 
elling of the spectral energy distributions, galaxy colours, 
discussion of the galaxy populations and source counts in 
the shorter Spitzer-IKAC bands and the longer wavelength 
MIPS bands (70/im & 160pm) will be given in a later contri- 
bution IPearson et al. 200511 . In section |21 two phenomeno- 
logical galaxy evolution models used to model the midin- 
frared counts are presented. In section |3 we discuss the re- 
sulting fits to the source counts at 15pm and 24pm. The 
conclusions are given in section |1] . A concordance cosmol- 
ogy of Ho = 72kms~^Mpc~'^,Q = 0.3, A = 0.7 is assumed. 



2 THE MID-INFRARED POPULATION 

The models presented are referred to as the Bright End evo- 
lution model and the Burst evolution model and are modi- 
fied and updated realizations of those presented in Pearson 
12001all . Pearson & Rowan- Robinson (IB^ . The previous 
models were based on the far-infrared 60pm luminosity func- 
tion of PSCz 77?^ 5 galaxies llSaunders et al. 2000l . The new 
models use the 15pm type dependent luminosity functions 
derived from the JS'O-ELAIS survey of Pozzi et al. lUTHUl 
to represent the normal -I- starburst/LIG/ULIG populations 
and Matute et al. 120041 to represent the Type-1 AGN (in- 
cluding Seyfert 1). Since Matute et al. 1200411 only provide 
a parametric luminosity function for the type 1 population, 
the 12pm luminosity function of Rush et al. (1993) is re- 
tained for the type 2 population. The use of mid-infrared 
luminosity functions as opposed to the more traditional far- 
infrared variants to model the source counts at 15 & 24pm 
minimize any errors due to assumed SED shape and K- 



corrections. The luminosity function of Pozzi et al. 12004t 
is used in preference to the 12pm & 25pm IRAS luminosity 
functions of Xu et al. 1200111 or Fan et al. 1199811 because 
the former is segregated specifically into normal quiescent 
and starburst populations based on ISO colour criteria . 

The adopted SED is the most important factor in the 
models. The complicated nature of the mid-infrared spectra 
and the variation between models, even for more normal 
galaxies can have a dramatic effect on the derived source 
counts of a galaxy population. In the final analysis of the 
ELAIS survey, Rowan-Robinson et al. 1200411 have shown 
that the infrared galaxy population can be divided into 4 
distinct spectral classes of normal quiescent (cirrus) galaxies, 
starbust (M82) galaxies, LIG/ULIG {IgLiR > 
sources (ARP220) and AGN (dust torus). 

The original normal galaxy used in Pearson 12001all has 
been replaced by the SED of Dale et al. 1200111 which incor- 
porates a more recent dust model (see Pearson & Efstathiou 
1200511 for a comaprison of cirrus galaxy SEDs). Although, 
for the normal galaxies, the effect of changing the SED in 
the 24pm band is minimal, at 15pm the original model pro- 
duced excess power in the differential counts at the level of a 
few mjy due to extremely strong and narrow PAH emission, 
masking the evolutionary upturn in the differential counts 
at that level. 

For the starburst template, the critical point in the SED 
is where the Wien-side of the dust emmission curve meets 
the mid-infrared PAH emmission. Analysis of model star- 
burst SEDs shows a large variation in resulting K-corrections 
due to 2 factors, the slope of the SED at this point (an in- 
creasing slope to shorter wavelengths producing a favourable 
negative K-correction) and the inclusion of distinct 12.7pm 
& 11.3pm features to produce a more accurate distribution 
of mid-infrared features (Boulanger et al. 1199811 . Dieden- 
hoven et al. (2001)). Thus we have retained the original 
Starburst (M82 like) SED of Efstathiou, Rowan- Robinson 
& Siebenmorgen 12000at but have replaced the SED in the 
range from 18-5pm, with the observed ISOCAM Circular 
Variable Filter spectrum of Forster Schreiber et al. 1200311 . 
Note that a similar adjustment was necessary in the models 
of Lagache et al. 1200411 in order to produce a reasonable 
fit to the Spitzer 24pm counts. 

The luminous infrared galaxy population is represented 
by the archetypical ARP220 SED of Efstathiou, Rowan- 
Robinson & Siebenmorgen 12000all which gives a good 
fit of the spectrum of ARP220 from sub-mm to optical 
wavelengths and is in reasonable agreement with the re- 
sults of Spoon et al. (i2004bl . Note that although it 
is known that ARP220 has an atypical far-infrared/mid- 
infrared luminosity ratio, we justify its use by the fact 
that it provides good fits to many of the ELAIS sources 
IRowan- Robinson 20041 and has more typical 6. 7- 12- 15pm 
luminosity ratios lElbaz et al. 2002ll . 

The AGN template uses the dust torus model of Rowan- 
Robinson 11995|l extended to near-infrared-optical wave- 
lengths by King & Rowan- Robinson 120031 . All the assumed 
SEDs are found to be consistent with the /SO-ELAIS colours 
of Rowan- Robinson et al. |2004j. For the calculation of the 
source counts the SEDs are smoothed by the respective fil- 
ter bandwidths corresponding to the /S'0-LW3 15pm and 
Spitzer-MIPS 24pm bands. 

We investigate two particular models, categorized by 
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their dominant populations of starburst (M82 like) and 
ULIG (Arp 220 like) sources respectively. The Bright End 
Model broadly follows the evolutionary scenario of Pearson 
& Rowan- Robinson (19961 assuming the starburst/ AGN 
and LIG/ULIG population luminosity function $(L, z) 
evolves in both luminosity, k, and density, g, as $(!//( 1 + 
z)'',0).{l + zy, with fc = 3.3 & 2.5 and 5 = 3.3 & 3.5 for the 
starburst and LIG/ULIG components respectively to z~l 
and constant thereafter. The Burst Model broadly follows 
the evolutionary scenario of Pearson 12U01at and assumes 
a similar power law evolution in luminosity and density, 
for the starburst and AGN populations and an initial ex- 
ponential burst -I- power law evolutionary scenario for the 
ULIG population. In this scenario, the starburst population 
evolves as (1 + z)'''^ , k = g — 3.2 while the more luminous 
galaxies evolve exponentially to z~l in density equivalent 
to (/ ~ 7 and then as a lower power law in both density 
and luminosity to higher redshifts. This violent evolutionary 
phase from z~0-l is consistent with extrapolations from the 
IRAS fednt source survey (Bertin et al. 119971 . Kim & Saun- 
ders 1199811 1. more extreme models of evolution in the ISO 
population | |Lagache, Dole, Puget 2003| l, models of bright 
sub-mm SCUBA galaxies (Small et al. 11997L Hughes et 
al. (troggt . Scott et al. Mortier et al. FMM ) and 

the most recent models of the Spitzer far-infrared poulation 
l |Lagache et al. 2004| |. In both models the normal galaxies 
are assumed to be non-evolving. 



3 COUNTS AND REDSHIFT DISTRIBUTIONS 

In figure the model fits to the observed counts in the 
ISO 15/xm and Spitzer 24/nm wavebands are shown for the 
Bright End (a) & Burst (b) evolutionary scenarios respec- 
tively. The counts are shown in Euclidean normalized differ- 
ential form. The two models are capable of simultaneously 
fitting the differential source counts in both the ISO Ibfim 
and Spitzer 24/im wavebands. The fits at 15/xm confirm the 
strong evolution to z~l implied by the ISO surveys while 
the counts at fainter fluxes at 24/im indicate a population 
of sources at higher redshift. The more conventional Bright 
End model predicts that star-forming galaxies rather than 
ULIG sources dominate the 15/im source counts. 

Pozzi et al. 12UU4II have fitted the ELAIS-S 15/im 
counts assuming an evloving starburst population with 
(1 -I- z)'°'®, k — 3.5, g — 3.8. Although assuming the Pozzi et 
al. 12UU4II evolution for the starburst component presented 
in this work does indeed provide a good fit at 15/im, it also 
produces an excess at around ImJy in the 24/im band, thus 
counts at fluxes of around 0.5-2mJy, in the Spitzer 24/im 
band can provide good constraints on the z<l evolution. 

In contrast, the more extreme burst model predicts that 
the upturn at mjy level fluxes at 15/im and the peak of the 
counts at 24/xm are attributed to the emergence of a new 
population of LIG/ULIG sources. Finally, both the Bright 
End Model and the Burst Model predict that the peak in 
the 24/im differential counts at ~0.2mJy is produced by the 
luminous LIG-ULIG population consistent with the early 
analysis of the Spitzer 24/im data (Yan et al. 12UU4II . Le 
Floc'h et al. TIUU^ ) . 

Thus, although from figure both models are in agree- 
ment as to the origin of the dominant source population at 



24/im, they difi'er in their prediction for the source of the up- 
turn (note, not the peak) in the differential counts at 15/tm, 
thus it is the counts at 15/im that have a better chance of dis- 
tinguishing between the two evolutionary scenarios. For the 
Burst model, we would expect to see a significantly higher 
fractional contribution of LIG/ULIG sources (an order of 
magnitude in the difi^erential counts at the mJy level) that 
simply are not appreciably present in the Bright End model. 

Note that at the depth of the ELAIS survey {S15 > 
Q.lmJy) the predictions of both models are consistent with 
the near-infrared and optical imaging of the ELAIS-S field 
by Pozzi et al. 12UU3II . with the majority of 15/im sources 
being luminous starburst galaxies and approximately 20% of 
sources attributed to AGN (Type 1 & 2). Rowan- Robinson 
et al. 1200411 find a high proportion ~14% of ULIGs in 
the ELIAS final band merged catalogue. This proportion 
is ~10% for the Bright End model. The Burst model pre- 
dicts a fraction of ~30% which is still consistent with the 
ELAIS results if the blank optical identifications (~10%) can 
be attributed to luminous infrared galaxies. Spectroscopy 
of the ELAIS-S sources has been made by La Franca et 
al. (2004 ). Of ~87% of sources spectroscopically identified, 
two populations dominated the 15/im sources over the range 
0.5-150mJy, with approximately three quarters contribution 
coming from normal to extreme starbursts and the remain- 
ing contribution coming from AGN. Around ~18% of the 
ELAIS-S sources have no counterpart to optical maginitudes 
R>23. Moreover La Franca et al. 120041 have shown that 
it is these sources are responsible for the upturn in the dif- 
ferential sources counts around ~lmjy at 15/im and that 
they are expected to have luminosities in the 10^^ — IO^^Lq 
range, i.e. luminous and ultra-luminous infra-red galaxies. 

From figure it can be seen that although both the 
fainter 15/im & 24/im differential counts can be well fitted 
at the fainter fluxes, the brighter counts are slightly un- 
derestimated in the Spitzer 24/im band. This could be at- 
tributed to some variance in the Spitzer fields since the initial 
counts from the SWIRE survey are somewhat lower at fluxes 
brighter than lOmJy (Shupe et al., in preparation). Alterna- 
tively, this effect could be due to a low normalization in the 
assumed ELAIS-S field luminosity function since although 
there is a good agreement with the evolution found in the 
deeper ISO surveys, including the upturn in the differential 
counts, in the brighter flux range of 0.5-lOmJy, the counts 
in the ELAIS-S fleld are systematically lower than the (cor- 
rected) Northern ELAIS flelds and the fainter ISO surveys 
(Serjeant et al. 120001 . Vaisanen et al. 120021 . Vaccari et 
al. (UnnUl, Rodighiero et al. i^OnHl 'l. 

Both evolutionary models are found to be broadly con- 
sistent with the 15/im redshift distribution at bright fluxes 
with a median redshift consistent with that found for the 
ELAIS sources (Pozzi et al. 120041 . Rowan- Robinson et al. 
(2004)). However, we predict an additional component at 
higher redshift ~1 in line with the assumptions of La Franca 
et al. 12004|l but not present in the estimations of Rowan- 
Robinson et al. (|2'04 ) . In flgure |5| the number-redshift 
distributions at fainter fluxes corresponding to the approx- 
imate peaks in the differential counts at Sis = OAmJy & 
S24 ~ 0.25mJy are shown respectively for the ISOCAM 
15/im & Spitzer 24/im bands for the Bright End (panel a) 
& Burst (panel b) evolutionary models. At the flux corre- 
sponding to the peak in the 15/im differential counts we find 
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Figure 1. Model fits to source counts in the ISOCAM 15/jm and 
Spitzer 24/im bands for tlie Bright End (a) & Burst (b) mod- 
els. Counts are presented in differential form, normalized to the 
non-evolving Euclidian universe. Total counts are shown with in- 
dividual components corresponding to Normal (dash), Starburst 
(dot), LIG/ULIG (dot-dash) sources and AGN (triple dot-dash). 
ISO 15/jm source counts from faint to bright fiuxes are from the 
ISO cluster lens survey - Altieri et al. (1999), Hubble Deep Field 
- Oliver et al. (1997), Aussel et al. (1999)), Lockman Hole - El- 
baz et al. (1999), ELAIS-S - Gruppioni et al (2002), ELAIS-N - 
Serjeant et al. (2000) and normalized IRAS counts from 12;im - 
Rush et al. 1993. The Spitzer-MlPS 24/xm source counts are from 
Papovich et al. (2004). 



good agreement with the fainter ISO surveys that predict a 
median redshift between z^O. 8-1.0 and domination by pop- 
ulations of luminous infrared sources (Elbaz et al. J^OO^ . 
Oliver et al. (I2UU2L Mann et al. Ill997t . At this redshift 
(z~0.95) the 7.7/im rest frame feature passes through the 
ISO 15/im band producing strong K-corrections, which com- 
bined with the strong evolution produces the redshift peak. 
The effect of the 9.7 /im silicate absorption feature is also 
prominent at z~0.5 and is consistent with the drop in den- 
sity of galaxies observed in the ELAIS-S field at this red- 
shift ILa Franca et al. 200411 . At 24/im a bimodal redshift 
distribution is seen at z~l & ~2.5 corresponding to the pas- 
sage of the 6.2-7.7/im and 11.3-12.7/im features respectively. 
A similar bimodal distribution has also been suggested by 
Chary et al. 12004t although the models presented here 
predict a greater contribution from higher luminosity in- 
frared sources. A significant deficit of sources is seen at 
which can be attributed to the effect of the silicate 
feature at 9.7/im. Although the two evolutionary scenarios 
predict a different origin for the dominant population at 
low redshifts, both models are unanimous in predicting a 
new population of infra-red luminous sources in the red- 
shift range, z=l-3. The SCUBA sub-mm population would 
seem likely candidates for this population of sources, indeed 
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Figure 2. Model Number-Redshift distributions for ISOCAM 
15^im and Spitzer 24/im bands for the Bright End (a) & Burst 
(b) models. The distribution are per sq.deg. with redshift bin 
size 5z = 0.1. Total numbers are shown with individual com- 
ponents corresponding to Normal, Starburst, LIG/ULIG sources 
and AGN. The distributions are for 5i5 = O.AmJy & 524 = 
0.25mjy and correspond to the approximate peaks in the differ- 
ential counts given in figure 1. 



observations in the MIPS 24/im band of SCUBA sources 
in the Lockman Hole detected almost all of the sources in 
the range z~ 0.5-3.5 l(E garni et al. 2004| l. Furthermore, the 
redshift distribution of the high-z Spitzer sources and their 
predicted luminosity range (LIG-ULIG) correlates well with 
the median redshift of 2.2-2.5 found for sub-mm sources 
(Chapman et al. ll^iMl . J^HUH . Frayer et al. |2004|). Lu- 
minous/ultraluminous galaxies are responsible for the peak 
at around 0.25mJy in the 24/im differential counts in both 
models. At this flux level both models predict of the or- 
der of more than 1000 such sources per square degree and 
typical sub-mm fluxes at 850/im of 3-5mJy over the red- 
shift range 1-2.5, in very good agreement with the observed 
submillimetre source density in this flux range observed by 
SCUBA lIBlain et al. 1999t . 



4 CONCLUSIONS 

The two evolutionary models presented are capable of si- 
multaneously fitting both the 15^m ISO and 24/im Spitzer 
normalized differential source counts. Careful consideration 
has been given to the assumed distribution of PAH features. 
The ISO counts are well fitted by starburst or LIG-ULIG 
sources strongly evolving to redshift ~1, depending on the 
model, while both models are unanimous in predicting a 
new population of luminous infrared galaxies at z~2-3 de- 
tected by the Spitzer observations at 24/im, that were previ- 
ously unexpected from fits to the ISO counts. These higher 
redshift sources are interpreted as being the infrared coun- 
terparts to the luminous sub-mm population. The number- 
redshift distributions at 15/^m & 24/im, while confirming the 
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strong evolution in the infrared population, show a bimodal 
distribution dominated by PAH emission features between 
3-13^m and the silcate absorption feature at 9.7/im. 

Although the original model of Pearson i200Ial pro- 
vided an excellent fit to the source counts at 15/im and other 
ISO wavelengths, like other pre-Spitzer models it fails to re- 
produce the bump in the Spitzer 24/im counts at fainter 
fluxes (see Papovich et al ("20041 for a comparision of pre- 
Spitzer models). In summary, the main difference in the as- 
sumed evolution in the new model formulation and that 
of the original models is an increase in the evolutionary 
strength in the l<z<2 region. Alternative post-Spiteer evo- 
lutionary models of Gruppioni et al. 12UU5II and Lagache 
et al. H20U4t can be broadly categorized by the Bright End 
evolution model and the Burst evolution model respectively. 
All of these models have required some adjustment of the 
mid-infrared Starburst spectral energy distribution in order 
to fit the observations. Both the models of Lagache et al. 
and the Burst evolution model are characterized by extreme 
evolution in the high luminosity infrared galaxy population 
(i.e. the LIG-ULIGs). These models predict that the domi- 
nant population at 15^m will be starburst galaxies while at 
24/xm they will be sources with infrared luminosities of the 
order of a few 10^^ — at redshifts l<z<3. In con- 

trast the model of Gruppioni et al. 1^005 ) and similarly the 
Bright End evolution model presented here, propose more 
traditional power law evolution and that the more numer- 
ous population at both 15 and 24^m will be the starburst 
galaxies, in fact the model of Gruppioni et al. f20051 as- 
sumes that the contribution of more extreme IR galaxies 
(ULIGs) at both 15 and 24^m is negligible. 

Deeper observations over larger areas in the 15/im band 
and photometric redshifts using mid-infrared multiband ob- 
servations of PAH emission and the silicate absorption fea- 
tures can break the degeneracy in the models by revealing 
the dominant population at fainter Ibfim fluxes (Starburst, 
LIGs, ULIGs). Spitzer has the ability for imaging in this 
wavelength range (16 & 22/xm) by using the peak-up cam- 
era of the infrared spectrograph (IRS, Houck et al. 12UU4II '). 
and Charmandaris et al. (2004) have discussed the potential 
for estimating redshifts and selecting galaxies at z>1.5 from 
the mid-infrared features by using the IRS. Takagi & Pear- 
son ()2()()5|l have shown that a combination of mid-infrared 
filters can detect luminous infrared galaxies at z~l-2 due 
to the strong silicate absorption expected in these sources 
JArmus et al. 2004t. 

The Spitzer peak up imager is however, limited by its 
small field of view (0.9'xl.3'). Deep, larger area observations 
in this wavelength range will be available from the Japanese 
infrared mission ASTRO-F due for launch in late 2005 (Mu- 
rakami (dnnU, Pearson et aL if^nni^ l. The ASTB.O-F iocal 
plane instruments, include an infrared camera (IRC) with 
a field of view of >10'x 10' (almost 100 times the area of 
Spitzer's peak up imager, allowing coverage of the deepest 
Spitzer fields in a single pointing), covering the near- mid- 
infrared in 9 bands from 2-24/.tm and covering the Spitzer 
wavelength desert in 4 extra bands at 9, 11, 15, 20^m, mak- 
ing it the ideal complement to Spitzer. Deep surveys are 
planned over sub-few square degree scales with sensitivi- 
ties down to 0.04mJy at 15/im ijPearson et al. 2001b^ . i.e. 
deep enough to directly probe significant samples of the 
sources responsible for the upturn in the 15/im counts. Tak- 



agi & Pearson 1200511 predict of the order of 1000 candidate 
LIG/ULIGs could be detected per square degree using the 
silicate-break method which will be a large enough sample to 
break the degeneracy in the models via specific colour cuts 
utilizing ASTRO-F's larger number of filters in the critical 
wavelength range that samples the PAH emission in galax- 
ies from 0<z<2.5 which is will be only partially possible and 
difficult & time consuming to achieve with Spitzer's peak up 
imager. 

In the next few years, multiwaveband observations with 
Spitzer and ASTRO-F pmhing different redshift regimes will 
examine the infrared active epoch Z'~0.5-1.5 (Franceschini et 
al. (2003), Elbaz et al. (120041 ') and the link between the 
infrared and sub-mm populations at z>2 IBlain et al. 20041 
providing a definitive record of the cosmic star formation 
history and cosmic fossil background of the Universe. 
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